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The  gas  turbine  combustion  system  design  and 
development  effort  Is  an  engineering  exercise  to  obtain 
an  acceptable  solution  to  the  conflicting  design  trade¬ 
offs  between;  combustion  efficiency,  gaseous  emis¬ 
sions,  smoke,  ignition,  restart,  lean  blowout,  burner 
exit  temperature  quality,  structural  durability,  and 
lire  cycle  cost.  For  many  years,  these  combustor 
design  trade-offs  have  been  carried  out  with  the  help 
of  fundamental  reasoning  and  extensive  component  and 
bench  testing,  backed  by  empirical  and  experience 
correlations . 

Recent  advances  in  the  capability  of  computational 
fluid  dynamics  (CFD)  codes  have  led  to  their  applica¬ 
tion  to  complex  three-dimensional  flows  such  as  those 
in  the  gas  turbine  combustor.  A  number  of  U.S.  Govern¬ 
ment  and  industry  sponsored  programs  have  made  signifi¬ 
cant  contributions  to  the  formulation,  development, 
and  verification  of  an  analytical  combustor  design 
methodologv  which  will  better  define  the  aerothermal 
loads  in  a  combustor,  and  be  a  valuable  tool  for  design 
of  future  combustion  systems.  The  contributions  made 
by  NASA  Hot  Section  Technology  (HOST)  sponsored  Aero¬ 
thermal  Modeling  and  supporting  programs  are  described 
in  this  paper. 


INTRODUCTION 


The  goal  of  gas  turbine  combustion  system  design 
and  deve’opment  is  to  obtain  an  acceptable  solution  to 
the  conflicting  design  trade-offs  between  combustion 
efficiency,  gaseous  emissions,  smoke,  ignition, 
restart,  lean  blowout,  burner  exit  temperature  quality, 
structural  durability,  and  life  cycle  cost.  For  many 
years,  these  combustor  design  trade-offs  have  been  car¬ 
ried  out  with  the  help  of  fundamental  reasoning  and 
extensive  component  and  bench  testing,  backed  by  empir¬ 
ical  and  experience  correlations.  The  ultimate  goal 
has  been  to  develop  a  reliable  combustor  design  system 
that  can  provide  quantitatively  accurate  predictions 
of  the  complex  combustion  flow  field  characteristics 
30  til'll  an  optimum  compustiOn  system  design 
can  be  achieved  within  reasonable  cost  and  schedule 
constraints . 
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Empirically  based  procedures  have  led  to  success¬ 
ful  evolutionary  combustor  improvements.  However,  as 
these  methods  are  experience-based,  they  are  not  well 
suited  when  combustor  design  requirement  are  signif¬ 
icantly  different  from  that  of  current  technology 
engines.  The  rapidly  developing  CFD  (Computational 
Fluid  Dynamics)  capability  is  providing  an  additional 
tool  in  the  design  process  which  can  have  a  powerful 
positive  influence  on  future  design  caoability.  In 
these  codes,  combustion  system  subcomponents  including 
diffusers,  fuel  injectors,  and  combustor  liners,  in 
addition  to  the  complex  internal  flow,  need  to  be  accu¬ 
rately  modelled.  To  achieve  this,  physical  sub-models 
and  accurate  numerical  schemes  must  be  developed  to 
describe  the  various  aerothermochemi ca 1  processes 
occurring  within  the  combustion  chamber. 

A  number  of  U.S.  Government  and  company  sponsored 
programs  have  made  significant  contributions  to  the 
formulation,  development,  and  verification  of  an 
analytical  combustor  design  methodology.  These  have 
included:  U.S.  Army  Combustor  Design  Criteria  Valida¬ 
tion  (Bruce  et  al.,  1979;  Mongia  et  al.,  1979,  Mongia 
and  Reynolds,  1979),  NASA  Swirling  Reci rculating  Flow 
(Srinivasan  and  Mongia,  1980),  NASA  Soot  and  NCx  Emis¬ 
sions  Prediction  (Snvatsa.  1980),  NASA  Primary  Zone 
Study  (Sullivan  et  al.,  1983),  NASA  Mass  and  Momentum 
Transfer  (Johnson  and  Bennett,  1981;  Roback  and 
Johnson.  1983;  Johnson  et  al . ,  1984),  NASA  Lateral  Jet 
Injection  (Lilley,  1986;  Ferrell  and  Lilley,  1985; 
McMurray  and  Lilley,  1986;  Ong  and  Lilley,  1986),  NASA 
Dilution  let  Mixing  (Srinivasan  et  al . ,  1982,  1984, 
1985;  Siinivasan  and  White,  1986,  Holdeman  et  al., 

1984;  Holdeman  and  Srinivasan,  1986;  Holdeman  et  al., 
1987a),  NASA  Transition  Mixing  Study  (Reynolds  and 
White,  1986;  Holdeman  et  al . ,  1987b),  NASA  HOST  Aero¬ 
thermal  Modeling  (<enworthy  et  al.,  1983;  Sturgess, 
1983;  Srinivasan  et  al . ,  1983a,  1983b),  NASA  Error 
Reduction  (Syed  et  al.,  1985),  industry  IR  &  D  pro¬ 
grams,  and  advanced  combustor  development  programs. 

The  NASA  Hot  Section  Technology  (HOST)  Combustion 
Program  has  supported  several  of  these  P’^ograms.  The 
overall  ohiecrivo  of  the  HPS'^  romhust-iAn  Projec*'  's 
develop  and  verify  advanced  analytical  methods  to 
iinpiove  the  capability  to  design  combustion  systems 
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for  advanced  aircraft  gas  turbine  engines.  This  objec¬ 
tive  is  being  approached  both  computationally  and 
experimental ly . 

Computationally,  HOST  first  sponsored  studies  to 
assess  and  evaluate  the  capabilities  of  existing 
aerothermal  models  (circa  1982).  Based  on  the  results 
of  these  assessments  and  other  studies  in  the  liter¬ 
ature,  HOST  supported  several  studies  to  develop  new 
and  improved  numerical  methods  for  the  analysis  of  tur¬ 
bulent  viscous  recirculating  flows,  with  emphasis  on 
accuracy  and  speed  of  solution. 

The  objectives  of  HOST  sponsored  experimental 
studies  were  to  improve  understanding  of  the  flow  phys¬ 
ics  and  chemistry  in  constituent  flows,  and  to  obtain 
fully-specified,  benchmark-qual ity  experimental  data 
suitable  for  the  assessment  of  the  capabilities  of 
advanced  computational  codes. 

This  paper  reviews  the  advances  in  the 
state-of-th3-art  in  combustor  aerothermal  modeling, 
while  highlighting  the  programs  supported  by  the  HOST 
Project  (Turbine  Engine  Hot  Section  Technology,  1982, 
1983,  1184.  1985.  1986,  1937).  Due  to  loagi.h  I'mita- 
tlons  not  all  programs  that  received  HOST  support  are 
included,  and.  for  completeness,  some  programs  that 
made  a  significant  contribution,  but  which  did  not 
draw  their  primary  support  from  HOST  are  discussed. 

AEROTHERMAL  MODELING  ASSESSMENT 

Gas  turbine  combustion  models  include  submodels 
of  turbulence,  chemical  kinetics,  turbulence/chemi stry 
interaction,  spray  dynamics,  evaporation/combustion, 
radiation,  and  soot  formation  and  oxidation.  A  very 
extensive  assessment  of  numerics,  physical  submodels, 
and  the  suitability  of  the  available  data  was  made  by 
three  contractors  under  Phase  1  of  the  HOST  Aerothermal 
Modeling  program  (Kenworthy  et  al.,  1983;  Sturgess, 
1983;  Srinivasan  et  ai.,  1983a,  1983b).  These  investi¬ 
gations  surveyed  and  assessed  current  models  and  Iden¬ 
tified  model  deficiencies  through  comparison  between 
calculated  and  measured  quantities.  Results  of  the 
assessment  by  Srinivasan  et  al . ,  (1983a,  1983b)  are 
summarized  by  Mongia  et  al .  (1986).  The  constituent 
flows  examined  included:  (1)  simple  flows  with  no 
streamline  curvature,  (2)  comple*  flows  without  swirl, 
and  (3)  complex  flows  with  swirl.  Geometries  for 
several  test  cases  from  each  of  these  categories  are 
shown  in  Fig.  2. 

k-e  Turbulence  Model 

The  k-e  model  is  the  simplest  turbulence  model 
that  is  suitable  for  recirculating  flow  calculations. 
This  model  achieves  closure  by  using  a  gradient  trans¬ 
port  model  for  Reynolds  stress  with  an  isotropic  eddy 
viscosity,  for  flows  where  the  Isotropic  eddy  vis¬ 
cosity  assumption  is  not  valid,  the  k-c  model  may  be 
either  modified  (e.g.  low  Reynolds  number  correction, 
Richardson  number  correction)  or  replaced  with  an  alge¬ 
braic  or  differential  Reynolds  stress  mode'. 

Assessment  of  the  k-c  imodeKs)  of  turbulence 
showed  that  these  models: 

(1)  require  low  Reynolds  Number  correction  for 
predicting  wall  shear  flows,  and  streamline 
curvature  modifications  for  accurately 
predicting  curved  boundary  layers 

(2)  give  quantitatively  good  correlation  with 
data  for  simple  flows  and  non-recirculating 
swirl ing  flows 

(3)  give  quanti tat ' ve I y  reasonaDle  results  for 
nonswirling  recirculating  flows 


(4)  give  quantitatively  unsatisfactory 
correlation  with  data  for  complex  swirling 
flows  with  recirculation  zones 

(5)  give  quantitatively  unsatisfactory 
correlation,  but  predict  trends  correctly,  for 
complex  three-dimensional  flows. 

Algebraic  Stress  Model  and  Its  Modifications 

Mean  flow  predictions  with  this  model  agreed  with 
the  data  as  well  as  the  k-c  model  results,  therefore 
the  conclusions  above  also  apply  to  this  model.  In 
addition,  the  Algebraic  Stress  Model  gives  reasonable 
predictions  for  the  Reynolds  stress  components,  con¬ 
sistent  with  the  strengths  and  limitations  of  the  k-e 
models  (Mongia  et  al.,  1986). 

The  results  of  standard  k-c  and  algebraic  and  dif¬ 
ferential  Reynolds  stress  turbulence  models,  have  been 
compared  in  several  continuing  assessment  studies.  An 
example  comparison  (Mongia,  1987)  of  data  and  calcula¬ 
tions  using  a  hybri d/SIMPLE  numerical  scheme  is  shown 
in  Fig.  3.  This  flow  is  that  of  co-annular  turbulent 
jets  flowing  into  an  axisymmetric  sudden  expansion 
(Roback  and  Johnson,  1983).  In  this  figure,  velocity 
profiles  are  shown  at  downstream,  distance  from  0.11 
to  2.5  pipe  diameters  from  the  expansion. 

Scalar  Transport  Model 

Mongia  et  al.,  (1986)  reported  that  the  k-c  model 
with  specified  Prandtl  number  predicts  scalar  fluxes 
reasonably  well  for  flow  where  the  gradient  diffusion 
approximation  is  valid.  An  alternative,  the  algebraic 
scalar  transport  model,  has  the  capability  to  improve 
predictions  over  the  k-c  approach,  but  further  work  is 
needed  to  establish  its  validity  for  swirling  recircu¬ 
lating  flows. 

Turbul ence/Chemi stry  Interaction  Models 

It  was  also  concluded  by  Mongia  et  al.,  (1986) 
that  both  2-  and  4-step  reaction  schemes  showed  prom¬ 
ise  for  application  in  gas  turbine  combustors,  but 
need  to  be  further  val'd?i'ed  against  (data  simple 
flames.  The  modified  eddy  breakup  model  predicted 
trends  well,  and  it  was  recommended  that  it  should  be 
pursued  because  this  approach  could  be  easily  extended 
to  multistep  kinetic  schemes. 

Numerical  Accuracy 

A  significant  deficiency  identified  in  the  assess¬ 
ments  was  that  for  many  flows  of  interest  the  accuracy 
of  the  calculation  was  limited  by  the  numerical  approx¬ 
imations,  wherein  the  false  diffusion  is  of  the  same 
order  of  magnitude  as  the  turbulent  diffusion.  This 
masked  the  differences  between  turbulence  models  such 
that  very  different  (models  gave  essentially  the  same 
result,  and  sometimes  resulted  In  undeservedly  good 
agreement  between  data  and  predictions. 

K  false  diffusion  is  present,  the  numeiical  solu¬ 
tion  obtained  for  any  given  flow  depends  on  the  grid 
density  and  distribution.  An  example  of  the  compari¬ 
sons  made  in  the  assessment  program  is  given  by  the 
comparison  in  Figs.  4  and  5  between  measured  and  calcu¬ 
lated  temperature  distributions  downstream  from  a  row 
of  jets  entering  a  confined  crossflow.  This  flow  is  a 
constituent  flow  in  most  gas  turbine  combustors,  and 
has  been  treated  extensively  in  the  literature,  includ¬ 
ing  the  recently  completed  NASA  Dilution  Jet  Mixing 
program,  from  which  data  were  compared  with  three- 
dimensional  calculations  in  the  Phase  I  assessment 
study  by  Srinivasan  et  al.,  (1983). 
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The  calculated  and  experimental  results  shown  are 
for  a  single  row  of  jets  with  an  orifice  spacing  to 
diameter  ratio,  S/0,  =»  2  Injected  Into  a  ducted  main¬ 
stream  with  a  duct  height  to  orifice  diameter  ratio 
H/0,  «  8.  The  jet-to-ma1 nstream  momentum  flux  ratio, 

J,  for  this  test  was  25.32.  Calculations  for  this 
case  made  with  45x26x17  (19890)  nodes,  are  shown  in 
Fig.  4.  The  parameter  plotted  In  these  figures  is  the 
dimensionless  mean  temperatui'e  difference  ratio,  THETA, 
where  THETA  *  (Tm  -  T)/(Tm  -  Tj).  The  predicted  jet 
penetration  and  mixing  are  less  than  that  shown  by  the 
data. 

The  calculation  shown  in  Fig.  4  used  49  nodes  to 
simulate  each  jet.  It  is  generally  not  possible  to  use 
this  many  grid  points  in  such  a  small  region;  as  few 
as  four  may  be  used  in  practice  for  each  jet.  To  simu¬ 
late  the  accuracy  of  this  approximation,  calculations 
were  performed  for  the  same  flow  and  geometric  condi¬ 
tions,  but  with  a  27x26x8  (5615)  grid.  These  coarse- 
grid  calculations  (Fig.  5)  are  in  much  better  agreement 
with  the  data  than  the  fine-grid  calculations.  These 
and  other  calculations  in  Srinivasan  et  al.,  (1983b) 
clearly  demonstrated  that  the  three-dimensional  calcu¬ 
lations  were  not  grid  independent. 

Conclusions  from  the  Assessments 

The  major  conclusion  in  the  HOST  Aerothermal 
Modeling  Phase  I  assessment  studies  by  Kenworthy 
et  al.  (1983),  Sturgess  (1983),  and  Srinivasan  et  al. 
(1983a,  1983b)  was  that  the  available  computational 
fluid  dynamics  (CFO)  codes  provided  a  useful  combustor 
design  tool.  Although  significant  advances  have  been 
made  in  the  development  and  validation  of  multidimen¬ 
sional  gas  turbine  combustion  calculation  procedures, 
the  codes  assessed  were  only  qualitatively  accurate, 
especially  for  complex  three-dimensional  flows,  and 
further  work  was  needed.  It  was  concluded  that  both  a 
significantly  improved  numerical  scheme  and  fully- 
specified  experimental  data  (i.e.  both  mean  and  turbu¬ 
lence  flowfield  quantities,  with  measured  boundary 
conditions)  for  complex  non-reacting  and  reacting 
constituent  flows  were  needed  before  various  emerging 
physical  sub-models  of  turbulence,  chemistry,  sprays, 
turbulence/chemistry  interactions,  soot  formation/ 
oxidation,  radiation,  and  heat  transfer  could  be  prop¬ 
erly  assessed. 

A  SECOND  GENERATION  MODEL 

The  first  generation  combustor  design  procedure 
outlined  by  Mongia  and  Smith  (1978)  has  been  very  use¬ 
ful  for  developing  several  combustors  (Mongia  et  al . , 
1986)  that  exhibited  significant  technology  advances. 
However,  in  addition  to  the  model  deficiencies  identi¬ 
fied  in  the  assessments,  there  were  several  parameters 
of  importance  in  gas  turbine  combustor  design  that  the 
analytical  models  could  not  predict;  e.g.  gaseous  emis¬ 
sions,  soot  format lO'i,  flame  blow-out  limits,  combus¬ 
tor  pattern  factor,  and  liner  heat  transfer  These 
parameters  were,  however,  successfully  predicted  by 
well-established  semi-analytical  correlations  developed 
by  Plee  and  Mellor  (1980),  Lebfevre  (1985),  and  their 
associates.  Therefore,  a  combustor  design  procedure 
that  could  be  applied  to  current  and  future  gas  turbine 
engines  was  implemented  that  makes  Mse  of  emp’''’cal 
design  concepts  and  employs  analytical  modeling  tools 
to  represent  various  combustion  processes  (Rizk  and 
Mongia,  1986;  Mongia,  1987). 

This  method  makes  use  of  multidimensional  models 
to  establish  liner  flowfield  features  and  combustion 
characteristics.  The  analytical  results  are  then 
integrated  with  semi-empirical  correlations  for 


performance  parameters  of  interest.  That  is,  flow 
field  and  geometric  parameters  that  are  needed  in  the 
empirical  equations,  such  as  combustion  volume  and  the 
fraction  of  air  participating  in  the  primary  combustion 
reaction,  are  provided  by  the  analytical  calculations. 

Satisfactory  agreement  with  experimental  data  has 
been  shown  (Rizk  and  Mongia,  1986)  for  emissions,  per¬ 
formance  and  heat  transfer.  The  combustor  for  which 
data  were  available,  and  for  which  calculations  were 
performed,  is  shown  schematically  in  Fig.  6.  A  typi¬ 
cal  comparison  between  data  and  predictions  for  CO, 
unburned  hydrocarbons,  NOx ,  soot  emissions,  combustion 
efficiency,  pattern  factor,  and  lean  blowout  are  shown 
in  Figs.  7(a)  to  (g)  respectively.  The  model  is  in 
good  agreement  with  the  data  over  the  enti'^e  see-levei 
engine  operating  range.  Calculated  liner  wall  tempeia- 
tures  for  both  the  inner  and  outer  walls  of  this  com¬ 
bustor  are  shown  in  Fig.  8  for  three  typical  z-planes 
along  k  *  5,  14,  and  23.  Here  k  denotes  nodal 
planes  along  the  combustor  circumferential  direction. 
Although  no  direct  comparison  with  liner  wall  tempera¬ 
ture  data  was  made,  the  predictions  look  reasonable. 

AEROTHERMAL  MODELING  PHASE  II 

Based  on  the  recommendations  of  the  Phase  I 
assessment  studies,  activities  in  Phase  II  of  the  HOST 
Aerothermal  Modeling  program  concentrated  on  developing 
improved  numerical  schemes,  and  collecting  completely- 
specified  data  for  nonreacting  single  and  two-pha'e 
swirling  and  nonswirling  flows.  The  programs  initiated 
were:  Improved  Numerical  Methods;  Flow  Interaction 
Experiment;  and  Fuel  Injector/Air  Swirl  Characteriza¬ 
tion.  The  first  of  these  is  a  prerequisite  to  further 
model  development,  and  the  data  obtained  in  the  latter 
two  studies  will  be  used  to  validate  advanced  models 
being  developed  independently. 

Improved  Numerical  Methods 

The  hybrid  finite  differencing  scheme  employed  in 
generally  available  combustor  codes  gives  excessive 
numerical  diffusion  errors  which  preclude  accurate 
quantitative  calculations.  In  response  to  this  defi¬ 
ciency,  HOST  supported  three  programs  with  the  primary 
objective  to  identify,  assess,  and  implement  improved 
solution  algorithms  applicable  to  analysis  of  turbu¬ 
lent  viscous  recirculating  flows.  Both  solution  accu¬ 
racy  and  solution  efficiency  were  addressed  (Turbine 
Engine  Hot  Section  Technology,  1985,  1986.  1987;  Turan 
and  VanDoormal ,  1987) , 

For  most  practical  problems,  a  central  diffe’'- 
encing  scheme  would  he  ideally  suited  if  it  were 
unconditionally  stable.  Central  ifferencing  is  a 
simple  second-order  scheme  which  ,5  easy  and  straight¬ 
forward  to  implement.  However,  for  grid  Peclet  num¬ 
bers  larger  than  2,  central  differencing  can  lead  to 
over-  and  under-shoots  and  is  unstable.  The  hybrid 
(central /upwi  ndx  scheme  ii  staole  f-/i  all  Peclet  num¬ 
bers,  but  suffers  from  excessive  false  diffusion.  An 
alternative  scheme,  named  CONDIF  (Controlled  Numerical 
Diffusion  with  Intern''’  Fec’hack)  (Runchal  et  al  .  , 

1986)  has  unconditionally  positive  coetticients  and 
still  maintains  the  essential  features  of  central  dif¬ 
ferencing  and  its  second-order  accuracy. 

CONDIF  uses  cen*-',;!  di  Tfe,  e.nc’ :',g  wl.C'  -t-  1 

Where  Pe  >  z  ano  the  dependent  variable  varies  monoton- 
ically,  a  modified  central  differencing  scheme  is 
used,  otherwise  upwind  differencing  is  used.  CONDIF 
employs  just  enough  numerical  diffusion  to  ensure  sta¬ 
bility  based  internally  on  the  field  distribution  of 


the  variable,  rather  than  switching  to  upwind  differ¬ 
encing  whenever  Pe  exceeds  2.  Since  upwinding  is  done 
at  relatively  few  grid  points,  CONOIF  essentially  main¬ 
tains  the  second-order  accuracy  of  central  differen¬ 
cing,  and  false  diffusion  is  substantially  reduced. 

Another  advanced  numerical  scheme,  called  flux- 
spline  (Patanlcar  et  a1.,  1987),  is  based  on  a  linear 
variation  of  total  flux  (convection  +  diffusion 
between  two  grid  points.  This  is  an  improvement  over 
the  assumption  of  uniform  flux  used  in  hybrid  schemes, 
and  leads  to  reduced  numerical  diffusion. 

Both  of  these  schemes  have  been  used  to  solve  a 
variety  of  analytical,  two-dimensional  laminar  and  tur¬ 
bulent  flows  (Runchal  et  al . ,  1  987;  Patankar  et  al . . 
1987).  As  an  example,  results  for  a  laminar  flow 
(Re  -  400)  in  a  square  driven  cavity  are  shown  in 
Fig.  9.  This  flow,  shown  schematically  in  part  a),  is 
characterized  by  a  strong  recirculation  zone  typical 
of  many  physical  situations.  The  problem  was  solved 
with  both  CONDIF  and  flux-spine  schemes  on  a  uniform 
22x22  grid  and  compared  with  the  exact  analytical  solu¬ 
tion  and  a  hybrid  solution  on  an  extremely  fine  82x82 
grid.  Velocity  profiles  at  the  midsection  of  the  cav¬ 
ity  are  shown  in  Fig.  9(b).  Both  advanced  schemes 
show  improvement  over  the  hybrid  calculation. 

An  attractive  feature  of  both  CONDIF  and  flux- 
spline  schemes  is  that  their  extension  to  three  dimen¬ 
sions  is  relatively  straight-forward.  The  resulting 
linear  differential  equations  involve  only  seven  points 
as  opposed  to  27  points  needed  in  many  skewed-upwind 
schemes  (Syed  et  al.,  1985). 

In  addition  to  the  need  for  improved  numerical 
accuracy,  there  is  a  need  for  improved  computational 
efficiency  for  a  given  level  of  accuracy.  Typically 
the  continuity  and  momentum  equations  are  solved  sepa¬ 
rately,  and  then  linked  through  iteration  of  the 
pressure  term;  e.g.  SIMPLE  (Semi-Implicit  Method  for 
Pressure  Linked  Equations).  "Modi? icationsT  such  as 
SIMPLER  and  PISOT  have  been  shown  to  Improve  computa¬ 
tional  efficiency.  Other  advanced  schemes  (Turbine 
Engine  Hot  Section  Technology,  1985,  1986,  1987;  Vanka, 
1987),  such  as  block  correction  techniques  and  direct 
solution  of  the  coupled  equations  have  been  proposed. 
Calculations  with  the  latter  coupled  with  the  flux- 
spine  technique  have  shown  a  speed  increase  by  a  fac¬ 
tor  of  15  for  a  calculation  of  turbulent  flow  over  a 
backward-facing  step  (Mongta,  1987). 

Gas  Phase  Experiments 

An  experimental  study  of  the  interactions  between 
the  combustor  and  diffuser  systems  (Srinivasan  and 
Thorp,  1987)  is  in  progress  to: 

(1)  Identify  the  mechanisms  and  magnitude  of 
aerodynamic  losses  in  various  sections  of  an 
annular  combustor-diffuser  system 

(2)  Determine  the  effects  of  geometric  changes  in 
the  prediffuser,  dome,  and  shroud  on  these 
losses 

(3)  Obtain  a  data  base  to  assess  current  and 
advanced  aerodynamic  computer  models  for 
predicting  these  complex  flowfields 

(4)  Upgrade  the  analytical  models  based  on  the 
experimental  data 

(5)  Design  and  test  advanced  diffuser  systems  to 
•orify  the  accuracy  of  the  upgraded  analj'lli.dl 
mode  1 

Another  study  in  progress  will  obtain  comprehen¬ 
sive  mean  and  turbulence  measurements  of  velocity  and 
species  concentration  in  a  three-dimensional  flow  model 
of  the  primary  zone  of  gas  turbine  combustion  chambers 
(Turbine  Engine  Hot  Section  Technology,  1985,  1986, 


1987).  The  flowfield  of  interest  is  the  interaction 
between  swirling  flow  and  lateral  jets  in  a  rectan¬ 
gular  channel  (Fig.  10).  The  mainstreams  flow  enters 
through  5  swirlers  with  the  transverse  jets  injected 
from  both  the  top  and  bottom  duct  walls  with  either  2 
or  4  jets  per  swirler  at  1/2  or  1  channel  height  down¬ 
stream  from  the  swirler. 

These  experiments  are  being  conducted  on  both  air 
and  water  mul t i pi e-swi r ler  rigs,  as  well  as  single 
swirler  and  swirling  jet  rigs.  Fifteen  cases  (combina¬ 
tions  of  swirl  and  jet  strength  and  location)  are  under 
test  using  laser  sheet  light  and  dye  water  flow  visual¬ 
ization,  and  detailed  velocity  and  scalar  mean  and  tur¬ 
bulence  LDV  measurements  are  being  made  in  the  air  rig.  * 

A  key  feature  of  this  program  is  comparison  of 
model  calculations  against  the  data  obtained  to  ensure 
that  the  data  are  complete  and  consistent,  and  satisfy 
the  boundary  condition  input  requirements  of  current 
three-dimensional  codes.  Calculations  were  performed 
using  a  three-dimensional  code  (Srivasta,  1980)  for  all 
test  cases  before  the  experiments  were  begun.  Data  and 
both  previous  and  advanced  model  calculations  are  being 
compared  as  data  are  obtained. 

Fue 1 -Injector /A  1 r-Swi r1  Characterization 

The  objective  of  this  study  is  to  obtain  fully- 
specified  mean  and  turbulence  measurements  of  both  gas 
and  droplet  phases  downstream  of  a  fuel  injector  and 
air  swirler  typical  of  those  used  in  gas  turbine  com¬ 
bustion  chambers. 

The  flowfield  of  interest  is  an  axisymmetric 
particle-laden  jet  flow  with  and  without  confinement 
and  co-annular  swirling  air  flow.  Approximately  30 
cases  are  under  test  with  both  glass-bead  particle¬ 
laden  jets  and  liquid  sprays,  with  various  combination 
of  swirl  strengths  and  confinement  (Turbine  Engine  Hot 
Section  Technology,  1985,  1986,  1987).  Measurements 
of  mean  and  turbulence  quantities,  for  both  gas  and 
solid  phases  are  being  made  using  a  2-component  Phase/ 

Doppler  LDV  particle  analyzer  (McOonell  et  al.,  1987). 

Calculations  were  performed  for  all  test  cases 
with  a  two-dimensional  TEACH-type  nonreacting  turbulent 
viscous  two-phase  flow  code  before  the  experiments  were 
begun.  Data  and  both  previous  and  advanced  model 
calculations  are  being  compared  as  data  are  obtained 
(Mostafa  et  al . ,  1987.  1988;  Nikjooy  et  al.,  1988). 

In  the  first  series  of  tests,  the  developing 
regions  of  unconfined  single  and  two-phase  flows,  with 
105  pm  glass  beads,  have  been  examined  experimentally 
and  analytically  for  particle-to-gas  mass  loadings  of 
0.2  and  1.0.  Data  and  calculations  for  the  latter  are 
shown  in  Fig.  11.  A  two-component  Phase/Doppler  sys¬ 
tem  was  used  to  map  the  flowfield,  including  particle 
number  density,  and  two  orthogonal  components  of  veloc¬ 
ity  for  both  phases. 

Calculations  are  shown  for  both  deterministic  and 
stochastic  treatments  of  the  particles,  using  a  two- 
phase  k-c  model.  Both  treatments  of  the  particles 
give  the  same  gas-phase  axial  velocity  profiles,  how¬ 
ever,  the  stochastic  approach,  which  attempts  to  model 
particle/gas  phase  interactions,  gives  better  agree¬ 
ment  for  particle  quantities  than  the  deterministic 
approach  which  ignores  turbulence  interactions. 

Another  experimental  program  was  conducted  to 
outaiii  i nfoi mat i or,  c''  the  characteristics  o'"  the  ipray 
produced  by  a  gas  turbine  fuel  injector  (McVey  et  al., 

1988a,  1988b).  The  objective  of  this  study  was  to 
obtain  spatially-resolved  information  on  both  the 
liquid  and  gaseous  phases  of  the  spray  flow  field  unUer 
conditions  of  high-flow,  high  velocity,  and  high  swirl 
that  are  typical  of  engine  operation.  Measurements 
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were  made  with  a  high-resolution  spray  patternator,  a 
two-component  laser  velocimeter,  and  a  single-component 
Phase/Ooppler  particle  analyzer. 

The  comprehensive  experimental  data  generated  in 
these  programs  will  be  used  to  validate  advanced  models 
of  turbulence,  scalar,  and  spray  transport,  including 
two-ecjuatton  turbulence  models,  algebraic  and  differen¬ 
tial  Reynolds  stress  models,  scalar  and  scalar-velocity 
transport  models,  and  Eulerian  and  Lagranglan  determin- 
*  istic  and  stochastic  spray  models. 

SUMMARY 

i  Although  significant  progress  has  been  made  in 

the  development  of  three-dimensional  analytical  CEO 
codes  and  their  application  in  future  gas  turbine  com¬ 
bustor  design,  these  codes  are  neither  sufficiently 
comprehensive  nor  quantitatively  accurate  enough  to 
permit  a  complete  design  alone.  They  are,  however,  a 
valuable  component  in  an  evolving  combustor  design 
methodology  in  which  their  capability  is  integrated 
with  the  substantial  base  of  empirical  experience  and 
one-dimensional  flow  model i ng . 

CONCLUDING  REMARKS 

The  NASA  HOST  sponsored  Aerothermal  Modeling 
Phase  II  programs  will  lead  to  significant  improve¬ 
ments  in  our  technical  ability  to  predict  nonreacting 
gas  turbine  combustor  flow  fields  with  and  without 
spray  injection.  Significantly  enhanced  capabilities 
for  accurately  predicting  combustor  aerothermal  per¬ 
formance  and  wail  temperature  levels  and  gradients 
will  require  further  improvements  in  numerical  schemes 
and  physical  submodels.  It  is  equally  important  to 
collect  fully-specified  reacting  flow  data,  similar  to 
what  is  being  done  for  nonreacting  flows  under  HOST 
Phase  II,  for  both  complex  constituent  flows,  and 
generic  gas  turbine  combustors. 

In  parallel,  work  should  continue  in  the  formula¬ 
tion  and  systematic  validation  of  turbulent  combustion 
models  for  reacting  sprays  and  multidimensional  heat 
transfer  models.  These  capabilities  will  provide  the 
tools  needed  to  analytically  conduct  the  combustion 
trade-off  studies  so  that  optimum  future  combustion 
systems  can  be  designed,  fabricated,  and  developed 
within  acceptable  cost  and  schedule  constraints. 
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SWIRLING  FLOW  IN  A  PIPE  EXPANSION 
JANJUA  ET  AL  (1982) 
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FIGURE  3,  -  COMPARISON  OF  MEASURED  MEAN  AXIAL  VELOCITY  PROFILES 
FOR  COANNULAR  JETS  DOWNSTREAM  OF  AN  AX  I  SYMMETRIC  SUDDEN  EX¬ 
PANSION,  WITH  Calculations  made  using  three  turbulence  models. 
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OF  COOL  JETS  INJECTED  INTO  A  CONSTANT-TEMPERATURE  CROSS  FLOW  FROM  THE  UPPER  WALL  OF  A  CONSTANT  AREA  DUCT  (J  =  25.32 
S/H  =0.25,  H/D  =  8,  19  890  NODES. 
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EIGUKE  S.  -  COMPARISON  BETWEEN  MEASURED  AND  CALCULATED  DIMENSIONLESS  TEMPERATURE  DIFFERENCE  RATIOS  DOWNSTREAM 
OF  COOL  JETS  INJECTED  INTO  A  CONSTANT-TEMPERATURE  CROSS  FLOW  FROM  THE  UPPER  WALL  OF  A  CONSTANT  AREA  DUCT  (J 
S/H  -  0.?5,  H/D  =  8,  5615  NODES). 


DISTANCE,  M 


FIGURE  8.  --  CALCULATED  LINER  WALL  TEM¬ 
PERATURES  AT  MAXIMUM  POWER  CONDITION 
FOR  COMBUSTOR  IN  FIG.  6. 


DRIVEN  CAVITY 


(B)  VELOCITY  PROFILES  AT  SECTION  A-A. 

FIGURE  9.  -  CALCULATIONS  OF  LAMINAR  FLOW  IN  A  SQUARE 
(2-D)  DRIVEN  CAVITY. 
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FLOW  CONTROL  VALVE 

FIGURE  10.  -  TEST  SECTION  GEOMETRY  FOR  EXPERIMENTAL  STUDY  OF  IN 
TERACTION  BETWEEN  FLOW  FROM  MULTIPLE  SWIRLERS  AND  TRANSVERSE 


(C)  PARTICLE  NUMBER  DENSITY. 

FIGURE  11.  -  RADIAL  PROFILES  OF  GAS-AND  SOLID-PHASE  MEAN 
FLOW  COMPONENTS  AND  PARTICLE  NUMBER  DENSITY  AT  A  PARTI¬ 
CLE  NUMBER  DENSITY  AT  A  PARTICLE-TO-GAS  MASS  LOADING 
RATIO  OF  1.0. 
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